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Abstract 

Using the ideas of supersymmetry and shape invariance we show that the eigenvalues 
and eigenfunctions of a wide class of noncentral potentials can be obtained in a closed 
form by the operator method. This generalization considerably extends the list of exactly 
solvable potentials for which the solution can be obtained algebraically in a simple and 
elegant manner. As an illustration, we discuss in detail the example of the potential 

.r^ n .^ 2 ^ C D F G 

V(r, 6, 6) = —r + ^+ , . + ^ + ^ . . o — r + 



4 r'^cos^O r'^si'n?9 si'nP'acf) r'^si'n?9 cos'^ac, 

with 7 parameters. Other algebraically solvable examples are also given. 
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I. Introduction 

Most text books on nonrelativistic quantum mechanics show how the harmonic 
oscillator problem may be elegantly solved by the raising and lowering operator method 
[1]. Sometime ago, using the ideas of supersymmetry [2] and shape invariant potentials 
(SIP) [3], a generalization of this method [4] was used to handle many more potentials 
of interest. It was found there that this enlarged list included essentially all the solvable 
potentials found in most text books on nonrelativistic quantum mechanics. It turns out 
that most of these potentials are either one dimensional or are central potentials which 
again are essentially one dimensional but on the half-line. It is natural to enquire if one 
can also solve some noncentral potentials in three dimensions by the operator method . 
The purpose of this paper is to show that a large number of noncentral potentials that are 
separable in their coordinates and may have upto seven parameters can be solved by the 
operator method by using the known results for the various SIP. More precisely, we show 
that the energy eigenvalues and eigenfunctions of noncentral but separable potentials in 
spherical coordinates r,$,(f) can be simply written down by applying supersymmetry and 
shape invariance successively to the (f), 9 and r dependent potentials. 

The procedure outlined here is elegant, and simple, so that a student of quantum 
mechanics should be able to appreciate it. Indeed, we feel that the material presented 
here, along with that given in ref [4] , may be profitably included in a graduate quantum 
mechanics course. Accordingly, we have kept this article at a pedagogical level and made 
it self-contained . In sec. II , we give a brief review of the main ideas of supersymmetric 
quantum mechanics and the concept of shape invariance, and then show how to obtain the 
energy eigenvalues and the eigenfunctions of such potentials by a generalized operator 
method. Sec. Ill is the core of this article. In it we explain how the results for the 
known SIP may be used to algebraically obtain in a closed form the eigenvalues and the 

2 



eigenfunctions for noncentral but separable potentials . As an illustration we discuss in 
detail the noncentral potential 



f n 2 ^ C D F G 

Vi(r,6',0) = — r + ^ + o ■ 2fl + ^ ^ + o ■ 2fl • 2 - + 9 ■ 2fl 



(1.1) 



wlierea;,^, C,D, G,F and a are arbitrary parameters, and show that the corresponding 
energy eigenvalues are 



n,ni,n2 



(2n2 + l) + (5 + /i)^ 



(1.2) 



where 



{2ni + l) + \lD + - + {C + 



2^ in2 



(1.3) 



The expressions for the energy eigenvalues and the eigenfunctions of the other noncentral 
potentials are also given. 

II. Supersymmetry and Shape Invariance 

Consider a potential V- (x) whose ground state wave function i/;q ^ = 'i/^q{x) is known 
and whose ground state energy has been adjusted (by adding a suitable constant to 
V-{x)) to -E'o"'' = 0. Then the Schrodinger equation for the ground state is (we shall use 
h = 2M = 1 throughout this article) 



H_iPo = i-^ + V-ix))iPo= 0, 



so that 



V_{x) = 



If we now define the operators 



A=^^+Wix),A- = -^^+Wix), 



(2.1) 



(2.2) 



(2.3) 



where the superpotential W(x) is related to ipo by 

Wix) = ; M^) = Nexp[- J Wiy)dy] , (2.4) 



then one finds that 



H_=A+A=-^ + W^-W'(x). (2.5) 



The supersymmetric partner potential is given by 

H+ = AA+ = --^ + W'^ + W'{x) . (2.6) 

If i/jn ^ and i/ji^^ denote the eigenfunctions of the Hamiltonian H_ and with eigen- 
value -En"^ and respectively then it is easily shown that 

Ei+^ = E^t , (2.7) 

V'i^) = , (2.8) 

V'i^i = (^i+^)-^/^A+V'i+) . (2.9) 

Further, the reflection and transmission coefficients (or phase shifts) of the two potentials 
are also related [2] . 

The underlying reason for this relationship is the fact that there is a supersymmetry 
in the problem. In particular it is easily shown that the supercharges Q and Q'^ defined 

along with the Hamiltonian H given by 
satisfy the algebra 

q2 = = Q+^ [H, Q] = = [H, Q+] (2.12) 
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QQ+ + Q+Q = H. (2.13) 

It must be emphasized here that the supersymmetry only gives the relationship 
between the eigenvalues and eigenfunctions of the two Hamiltonians but does not yield 
the actual spectrum. For that purpose one needs an extra integrability condition called 
the shape invariance condition which was introduced by Gendenshtein [3] . Stated simply, 
if the pair of supersymmetric partner potentials V± (x) defined above are similar in shape 
and differ only in the parameters which appear in them, then they are said to be shape 
invariant. More specifically, if V±{x; gq) satisfy the requirement 

V+{x;ao) = V_{x;ai)+R{ai), (2.14) 

where qq is a set of parameters and ai is an arbitrary function of ao{ai = /(oo)) and 
the reminder R{ai) is independent of x, then V±{x;aQ) are said to be shape invariant. 
In such a case the energy eigenvalue spectrum of the Hamiltonian is given by [3] 

n 

M = Yl ^M; 4"^ (ao) = , (2.15) 
fc=i 

with a/s = f^{ao), i.e., the function f applied k times. Subsequently, Dutt et al [5] showed 
that the eigenfunctions i/^n"^ of can also be written down algebraically 

ip^~\x;ao) = A+{x;ao)A+{x;ai)...A+{x;an-i)ipo^\x;an) (2.16) 

which is clearly a generalization of the operator method for the celebrated harmonic 
oscillator potential. Later on Dabrowska et al [6] algebraically obtained the explicit 
expressions for ip^^ for all the known SIP by using Eq.(2.16). 

It turned out [7] that there are twelve SIP for which ai = ao + a;a being a constant. 
Out of these , eleven are contained in the table in ref.[4], while the remaining one was 
pointed out by Levai [8] .Recently it has been realized [9] that out of these twelve, two are 
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not really independent, so that one really has ten independent SIP (with ai = qq + a), 
and for all of which and ipn can be obtained algebraically. Most of these potentials are 
also contained in the factorization approach of Schrodinger [10] . It is worth repeating that 
all the SIP are either one dimensional or central potentials which are again essentially 
one dimensional in nature. 

III. Shape Invariance and Noncentral Potentials 

Noncentral potentials are normally not discussed in most text books on quan- 
tum mechanics. This is presumably because most of them are not analytically solv- 
able. However, it is worth noting that there is a class of noncentral potentials in 3- 
dimensions [11] for which the Schrodinger equation is separable. In spherical coordinates 
r, ^, </) (0 < r < oo, < ^ < TT, < </) < 27r) , the most general potential for which the 
Schrodinger equation is separable is given by 



K(,,,) = VW + I<P + ^, (3.1) 



where V{r), V{9) and V{(f)) are arbitrary functions of their argument. Although we 
use the same symbol V for simplicity of notation,the functions V{r),V{9) , and F(</)) 
need not be the same. The crucial point is that for each of V{r),V{9) and V{4>) one 
can choose a well known SIP and hence can immediately write down the corresponding 
energy eigenvalues and eigenfunctions algebraically. Thus by applying supersymmetry 
and shape invariance successively, i.e. first to V{(j)), then to V{9) and finally to F(r), 
we can algebraically solve a large class of noncentral potentials. First , let us see why the 
Schrodinger equation with a potntial of the form given by Eq.(3.1) is separable in the 
(r, 9, (j)) coordinates. The equation for the wave function \E'(r, 9, (p) is 



Q2^ ' r dr' r^^d9^ 89' r^sin^^ dcjy' 
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{E-V)-^. (3.2) 



It is convenient to write ^{r, $, 0) as 

Substituting Eq.(3.3) in Eq.(3.2), we obtain 

1 d^i? 1 11 £h 1 2/,! 1 r 1 ^ 



i? dr2 ' ^ ^ 4:r^ r^^ H d9^ ^ ^ 4 ^ r2sin26' ^ K d(/>2 

(3.4) 



Suppose K{(f)) obeys the equation 

cPK 



+ V{4))K{(t)) = m^K{(t)) . (3.5) 



d02 

Substituting this in Eq.(3.4), the 0-dependence in it is eliminated, resulting in the equa- 
tion, 

^ "^'^ + [V{r) -^] + + V{e) + (m^ - i)cosec2^] = E . (3.6) 



Next, let iy(^) obey the equation 

A further substitution of this equation in Eq.(3.6) then completes the separation of the 
variables , giving the radial equation , 

j2 d (j2 _ 1\ 

+ [Vir) + R{r) = ER{r). (3.8) 

The three solutions given by Eqs. (3.5), (3.7), and (3.8) may be implemented algebraically 
by using the well-known results for the SIP in each case. 

As one example, we discuss the potential given by Eq.(l.l) in some detail. Compar- 
ing Eqs. (1.1) and (3.1) , we obtain 

Vi (0) = F cosec^acj) + G sec^ acj) , (3.9) 
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where the subscript 1 refers to the particular example under cosideration. Now, we know 
from [4] that if the superpotential W is chosen as 

W ^ Atanacf)- Bcotacf) , (3.10) 

then, using Eqs.(2.5) and (2.6) , the corresponding (shape invariant) potentials are 
given by 

F^(0) = -{A + Bf + A{A ^ a) sec^acf) + B{B ^ a) cosec^a^ . (3.11) 

Hence, 

V+{A, B, a, x) = V_{A + a,B + a,a,x) + {A + B + 2af -{A + Bf . (3.12) 

On comparing Eqs.(3.9) and (3.11) we see that , discounting the overall constant 
-{A + By, Vi and V_are identical if 

F = B{B-a), G = A{A-a). 

Thus the energy eigenvalues m\ of Eq.(3.5) with F((/>) = Vi((/)) follow from Eqs.(2.14), 
(2.15) and (3.12) , 

ml = {A + B + 2naf , n = 0, 1,2, ... (3.13) 
The corresponding eigenfunctions are given by [4] , 

Ki(0) = (sina0)^/" (cosq;(^)^/"P#"^'"~^ (cos2q;0) . (3.14) 

Having solved the eigenvalue equation for the 0— part, we proceed to solve the 
Schrodinger equation given by Eq.(3.7) for the ^-vavariable .Again , comparing Eqs.(l.l) 
and (3.1), we see that 

Vx{d) = Ccosec^e + Dsec^e . (3.15) 
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The resulting Schrodinger equation for Hi {$) is 



+ [ (C + mf - 1/4) cosec^^ + i:»sec^^] = i^i . (3.16) 



Using the same algebraic procedure as before , the eigenvalues if and the eigenfunctions 
Hi{e) of Eq.(3.16) are found to be 

ll = {Ai + Bi + 2ni + ^f, (3.17) 

Hi{d) = (sin^)^i+V2(^Qg^)Ai p^^/'^i-^(cos2^) , (3.18) 

where 

D = Ai(Ai-l), C + ml = Bl. (3.19) 

Notice that in general there is no degeneracy in li except in the special case when C = 
, and q: is a rational number. In particular for C = , is given by 



ll = l + 2nl + {2n + l)a+^/DTl/4:+^/FT^^+^/GT^^. (3.20) 

So if a = p/g , where p, q are integers with no common factors , then there is degeneracy. 
Finally, the radial Eq.(3.8) for the potential (1.1) is given by 



where E^^^ denotes the energy eigenvalue . Choosing the superpotential 



(3.21) 



W = A,r- , (3.22) 



the corresponding shape invariant potential V- is given by 



y_(r) = Alr^ + ^(^2 + 1) - ^2(2^2 + 3) . (3.23) 
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A little algebra will show that 



V+{A2,B2,r) = V_{A2,B2 + l,r)-A2{2B2 + 3)+A2{2B2 + l) . 

On comparing Eqs.(3.21) and (3.23) , we see that , apart from the overall constant 
—A2{2B2 + 3) , the radial potential in Eq.(3.21) is the same as V- if 

^ = Al 6 + ll = {B2 + l/2f . (3.24) 

The energy eigenvalues of the noncentral potential (1.1) may again be immediately writ- 
ten down as before . We find 

^il,n. = (2^2 + ^2 + ^) a;, (3.25) 
while the corresponding radial part of the eigenfunction is given by 

R,{r) = r^^+'e-^"^"/^ Ln!^^^ur^) . (3.26) 

Thus the total eigenfunction for the noncentral potential (1.1) is given by the form (3.3), 
with Ri{r), Hi{e) and Ki{(f)) given by Eqs.(3.26) , (3.18) and (3.14) respectively. 

In terms of the seven parameters, the energy eigenvalues are given by Eq.(1.2) 
(n, ni, n2 — 0, 1, 2, ...) , where we see the highly nontrivial dependence of E^^'^ on these 
parameters. Note that there is a degeneracy in the problem when either S and/or C are 
zero and a is a rational number. The maximum degeneracy is obtained in case both S 
and C are zero and a is a rational number and then the energy eigenvalues are given by 

EL%,n, = |2(ni+n2+na + l + |) + + 1/4 + ^^ + ay^+^/F + ayA^u; . (3.27) 

Since the degeneracy is usually associated with some symmetry in the problem, it would 
be interesting to enquire about the extra symmetry in the problem. For the case a = 1 
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(apart from 5 = C = 0), the symmetry is clear. In such a case , the potential is simply 
the spherically symmetric harmonic oscillator , with the perturbations of the form G/x^, 
Fjy^ and C j . However, in general, for any rational a. , the symmetry is not obvious. 

We thus have shown that by choosing F(r), Vid') and appropriately , the energy 
eigenvalues and the eigenfunctions of the noncentral potential (1.1) can be written down 
algebraically by using the results for the well known SIP. By going through the list of the 
well known SIP , one can now immediately write down E and for a class of noncentral 
but separable potentials. In particular notice that another possible form for the shape 
invariant V^(0) is [12] given by 

Vi (0) = Fcosec^CK^ + Gcota0 . (3.28) 

We may now use the relation [8] that if 

= - Acota0 - — , (3.29) 

then the corresponding (SIP) V- is given by 

V- = A{A - a) cosec^ + 2B cot - ^ . (3.30) 

Hence the energy eigenvalues are 

= (A + naf - f , (3.31) 
[A + na)'^ 

where 

A = {a+^/a'^ + AF)/2 ,B = G/2. (3.32) 

Note that in unbroken supersymmetry, A > , and the negative root is rejected.The 
corresponding eigenfunctions are found to be [8] 

K2{(l)) = (sina0)("+")e«"^P^-"-"-*'^'-"-"+*'^) (zcota0) , (3.33) 
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where 

B A 

s = A/a , X = — , a = —. r . (3.34) 

a{s + n) 

Similarly, the other possible forms for the shape invariant ^-dependent potential 

V{e) are 

V2 {9) = C cosec26' + D coiO , (3.35) 

and 

Vz{e) = Ccosec^e + i:>cosec^cot^ . (3.36) 

The energy eigenvalues and eigenfunctions of the Schrodinger eq.(3.11) for the potential 
V2{e) are [8] 

ll = {A, + m + - _f\i,2 ^ A'i = C + m\ (3.37) 

H2{e) = (sin^)(^^+^^+^)e«^^Pi:^^-"^-^-''^^'-^^-"^-^+''^^^(zcot^) , (3.38) 
where m is either m\ or 777,2 as given by eq.(3.13) or (3.31) respectively , and 



D 

2(^1 +^1 + I) 



ai = — ■ ^ . (3.39) 



In order to find the eigenvalues and the eigenfunctions for the potential V^{9), note 
that if 

W = -(Ai + ^)cot^-5icosec^ , (3.40) 
then the corresponding (SIP) V_(^) is given by 

V_{e) = {Al + Bl-h cosec^^ + 2Ai B^cotd cosec^ - (^1 + ]-f . (3.41) 

Hence the the eigenvalues and the eigenfunctions of the Schrodinger eq.(3.11) for the 
potential Vz{9) are [4,6], 

ll = {A^+n^ + \f , (3.42) 
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where 

+ C = Al + Bl,D^2AiBi , (3.44) 

with m being either mi or m2 as given by Eqs.(3.13) or (3.31) respectively. 

FinaUy, the other possible form for the shape invariant r-dependent potential V(r) 

is 

V2(r) = -- + A. (3.45) 



Observe that [4] if 



then the corresponding (SIP) V-{r) is given by 

l._M = _£! + ^H(^ + _£L_. (3.47) 

^ ^ r 4(^2 + 1)2 

Hence the energy eigenvalues and the eigenfunctions for the potential which follow from 
ref.[4,6] are 

4 

4(n2 + 52 + 1)2 ' ^^-^^^ 
i22(r) = y''^+'exp(-^y)L2f +i(y) , (3.49) 

where 

with Z being li or Z2 or Z3 depending on the form of V{9). We thus see that there are 
2 (/)-dependent, 3 ^-dependent and 2 r-dependent SIP. As a result, we can immediately 
construct twelve different noncentral potentials by taking various combinations of F(r), 
V{d) and V{(j)) and in each case we can write down the eigenvalues and the eigenfunctions 

13 



algebraically by using the well known results for the SIP as mentioned above. The main 
results can be summarised as follows 

(i) For six noncentral potentials, the energy eigenvalue are given by the formula (3.25) 
and for the remaining six cases by the formula (3.48) depending on if one is con- 
sidering the oscillator or the Coulomb potential. Ofcourse,the value of B2 will vary 
from potential to potential depending on the form of ^(6*) and V{(j)) as explained 
above. 

(ii) For all the 12 potentials, the eigenfunctions are given by 

tl,{T,eA)-Ri{r)H,{e)K,{ct>) 

where i,q = 1,2 and j = 1,2,3, depending on the choice of the potentials. For example, 
for the potential (1.1), V' = Ri{f)Hi{9)Ki{(f))^ where Hi and Ki are as given by 
eqs.(3.26), (3.18) and (3.14) respectively. 

(iii) In case V{6) is Vi{6) or V^{6) , there is a degeneracy in the spectrum for 5 = 
irrespective of the form of V{r) and V{(j)). On the other hand, when V{(j)) = Vi(0) 
there is a degeneracy in the spectrum if C = and a is a rational number .However, if 
V{(j)) = V2{(j)) and V{9) = ^2(6*) there is no degeneracy in the spectrum irrespective 
of the form of l/(r). Finally, in case ^(6*) = V^i(^) or ¥3(9) , and V{(j)) = Vi{(j)) , there 
is maximum degeneracy when both C and S are zero and a is a rational number for 
either form of V{r). Degeneracy in the spectrum is usually associated with some 
symmetry in the system , and it may be worthwhile to explore this connection. 

(iv) For the six noncentral potentials where V(r) is taken to be the Coulomb potential, 
one also has a continuous spectrum over and above the discrete one . It would be 
interesting to obtain the phase shifts for these noncentral potentials algebraically. 

IV. Conclusions 
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In this paper, we have shown that there exist twelve noncentral but separable poten- 
tials each of which have seven parameters and for each of which the eigenvalues and the 
eigenfunctions may be written down in a closed form algebraically using the well known 
results for the shape invariant potentials. We find this remarkable . In this context, it is 
worth remembering that even for central potentials, the most general solvable potentials 
are those of Natanzon type with six parameters but where the potential is not explicitly 
known [13]. 

Generalization of our technique to other noncentral potentials is quite straight for- 
ward. For example, one could instead consider the Schrodinger equation in the cylindrical 
coordinates and obtain other examples of the noncentral potentials with seven parame- 
ters for which the spectra can be written down by using the well known results for the 
SIP V{4)), V{p) and V{z), where p = \J -\- ip-. Other system of coordinates , or even 
higher dimensions may similarly be considered. 
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